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<-h Abstract 

The binding energies of a range of nuclei and hypernuclei with atomic number A < 4 and 
strangeness |s| < 2, including the deuteron, di-neutron, H-dibaryon, 3 He, ^He, 4 He, ^He, and 
A \He, are calculated in the limit of flavor-SU(3) symmetry at the physical strange-quark mass 

CN) with quantum chromodynamics (without electromagnetic interactions). The nuclear states are 

extracted from Lattice QCD calculations performed with n/ = 3 dynamical light quarks using 

T — I an isotropic clover discretization of the quark action in three lattice volumes of spatial extent 

L ~ 3.4 fm, 4.5 fm and 6.7 fm, and with a single lattice spacing b ~ 0.145 fm. 
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I. INTRODUCTION 



The structure and interactions of the light nuclei have been the focus of experimental and 
theoretical explorations since the infancy of nuclear physics. Yet more than one hundred 
years later, and despite having made remarkable progress in describing these systems in 
terms of nuclear forces that are well-constrained by experiment, we remain unable to predict 
the binding and interactions of any given nucleus with reliable estimates of the associated 
uncertainties. It has long been accepted that quantum chromodynamics (QCD) and the 
electroweak interactions produce the nuclear forces, and consequently are responsible for 
the structure and interactions of all nuclei. Unfortunately, the complexity of the QCD vac- 
uum has so far prevented the calculation of low-energy and medium-energy nuclear systems 
directly from QCD. Beyond recovering the results of decades of experimental investigation, 
it is crucial to establish and verify tools with which to perform such calculations, with quan- 
tifiable uncertainties, in order to determine the properties and structure of exotic nuclei, 
and of matter in extreme environments or in kinematic regimes where experiments are not 
possible or practical. 

The only known way with which to calculate the low-energy properties of hadronic and 
nuclear systems rigorously is Lattice QCD (LQCD). In LQCD calculations, the quark and 
gluon fields are defined on a discretized space-time of finite volume, and the path integral 
over the fields is evaluated numerically. While LQCD calculations deviate from those of QCD 
due to the finite distance between points of the grid (lattice spacing) and the finite volume 
of the grid (lattice volume), such deviations can be systematically removed by reducing 
the lattice spacing, increasing the lattice volume and extrapolating to the continuum and 
infinite-volume limits using the known dependences determined with effective field theory 
(EFT). Calculation of important quantities in nuclear physics using LQCD is only now 
becoming practical, with first calculations of simple multibaryon interactions being recently 
performed, although not at the physical values of the light-quark masses. Early exploratory 
quenched calculations of the nucleon-nucleon (NN) scattering lengths [1, 2] performed more 
than a decade ago have been superseded by n/ = 2 + 1 calculations within the last few 
years [3, 4] (and added to by further quenched calculations [5-7]). Further, the first quenched 
calculations of the deuteron [8], 3 He and 4 He [9] have been performed, along with rif = 2 + 1 
calculations of 3 He and multibaryon systems containing strange quarks [10]. In addition, 
efforts to explore nuclei and nuclear matter using the strong coupling limit of QCD have led 
to some interesting observations [11]. Recently, rif — 2 + 1 calculations [12, 13] and rif — 3 
calculations [14-16] have provided evidence that the H-dibaryon (with the quantum numbers 
of AA) is bound at a pion mass of ~ 390 MeV at the physical value of the strange-quark 
mass [12, 13] and over a range of SU(3)-degenerate light-quark masses with m n ~ 469 MeV 
to 1171 MeV [14, 15]. Extrapolations to the physical light-quark masses suggest that a 
weakly bound H-dibaryon or a near-threshold resonance exists in this channel [13, 17-19]. 
We have searched for bound states in other channels at ~ 390 MeV [12, 13], and evidence 
has been found for a bound state in the strangeness-4 system. This is consistent with 
model-dependent and EFT predictions of a bound state at the physical pion mass [20-22]. 
In addition to the identification of bound states, calculations of hyperon-nucleon scattering 
extrapolated to the physical pion mass (using leading-order (LO) EFT) have been performed 
and directly compared with the results of phase-shift analyses of experimental data [23]. 

In this work we focus on the lightest nuclei and hypernuclei and present results of the first 
LQCD calculations of a number of s-shell nuclei and hypernuclei with A < 5, including 3 He, 
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^He, iHe, ^He, and a five-body state AE°pnn in the limit of exact SU(3)-flavor symme- 
try (and consequently, exact isospin symmetry). Hypernuclear spectroscopy is enjoying an 
experimental renaissance with ongoing and planned programs at DA$NE, FAIR, Jefferson 
Lab, J-PARC and Mainz providing motivation for enhanced theoretical efforts (for a re- 
cent review, see Ref. [24]). Our LQCD calculations are performed using an isotropic clover 
quark action at the SU(3)-flavor symmetric point corresponding to the physical strange- 
quark mass, with m n = rriK = m v ~ 800 MeV. Three lattice volumes have been employed 
with spatial extent of L ~ 3.4 fm, 4.5 fm and 6.7 fm, and calculations of systems with 
non-zero total momentum (boosted systems) have been performed to investigate the volume 
dependence of binding energies [25, 26]. As this is the first calculation of hypernuclei with 
baryon number A > 2, it is prudent to establish benchmarks for future works. The spectra 
of nuclei will have the simplest structure at the SU(3) symmetry point, where the up, down 
and strange quarks have the same mass, allowing for a relatively uncomplicated analysis. 
While any common light-quark mass could have been used, the physical value of the strange 
quark mass was chosen so that only the (common) up and down quark masses deviated 
from their physical values, and also so that the four- and five-baryon systems would be 
well contained within the three selected lattice volumes. Further, such a large value of the 
pion mass, combined with the temporal extent of the gauge-field configurations, strongly 
suppresses thermal effects which are present in all calculations and can provide a systematic 
uncertainty in extracting the small energy differences present in nuclei. Only one rela- 
tively coarse lattice spacing, b ~ 0.145 fm has been used in the calculations, dictated by the 
available computational resources, and therefore, an extrapolation to the continuum has not 
been performed. Further, extrapolations to the physical light-quarks masses have not been 
attempted because the quark-mass dependences of the energy levels in the light nuclei are 
not known. Future calculations at smaller lattice spacings and at lighter quark masses will 
facilitate such extrapolations and lead to first predictions for the spectrum of light nuclei 
with completely quantified uncertainties that can be compared with experiment. 



II. LATTICE QCD CALCULATIONS 



A. Computational Overview 

Three ensembles of isotropic gauge-field configurations, generated with a tadpole-improved 
Luscher-Weisz gauge action and a clover fermion action [27], are used in this work. This 
particular lattice-action setup follows closely the anisotropic clover action of the ensembles 
generated by the JLab group that we have used in our previous calculations [4, 10, 12, 
13, 17, 23]. The parameter tuning and scaling properties of this action will be discussed 
elsewhere [30]. One level of stout smearing [31] with p = 0.125 and tadpole-improved tree- 
level clover coefficient c S w = 1-2493 are used in the gauge-field generation. Studies [28-30] 
of the PCAC relation in the Schrodinger functional indicate that this choice is consistent 
with vanishing 0(b) violations, leading to discretization effects that are essentially 0(b 2 ). 
The parameters of the ensembles are listed in Table I, and further details will be presented 
elsewhere [30]. As multibaryon systems are the focus of this work, relatively large lattice 
volumes are employed for the calculations, with correspondingly large values of m n L and 
m v T . In order to convert the calculated (binding) energies from lattice units (l.u.) into 
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TABLE I: Parameters of the ensembles of gauge-field configurations and of the measurements used 
in this work. The lattices have dimension L 3 x T, a lattice spacing b, and a bare quark mass b m q (in 
lattice units) generating a pion of mass m^. N SIC light-quark sources are used (as described in the 
text) to perform measurements on iV c fg configurations in each ensemble. The three uncertainties 
associated with the pion mass are statistical, fitting systematic and that associated with the lattice 
spacing, respectively. 



Label 


L/b T/b 


P 


b m q 


b [fm] 


L [fm] 


T [fm] 


77V [MeV] 


m-n-L 




iV cfg 


iV src 


A 


24 48 


6.1 


-0.2450 


0.145 


3.4 


6.7 


806.5(0.3)(0)(8.9) 


14.3 


28.5 


3822 


72 


B 


32 48 


6.1 


-0.2450 


0.145 


4.5 


6.7 


806.9(0.3)(0.5)(8.9) 


19.0 


28.5 


3050 


48 


C 


48 64 


6.1 


-0.2450 


0.145 


6.7 


9.0 


806.7(0.3)(0)(8.9) 


28.5 


38.0 


1905 


54 



physical units (MeV), a lattice spacing of b = 0.1453(16) fm has been determined for these 
ensembles of gauge-field configurations from the T spectrum [32] . 

The iV c fg gauge configurations in each of the ensembles are separated by at least 10 
HMC evolution trajectories to reduce autocorrelations, and an average of iV src measure- 
ments are performed on each configuration. The quark propagators were constructed with 
gauge-invariant Gaussian-smeared sources with stout-smeared gauge links. These sources 
are distributed over a grid, the center of which is randomly distributed within the lattice 
volume on each configuration, and the quark propagators are computed using the BiCGstab 
algorithm with a tolerance of 10~ 12 in double precision. The quark propagators, either un- 
smeared or smeared at the sink using the same parameters as used at the source, give rise 
to two sets of correlation functions for each combination of source and sink interpolating 
fields, labeled as SP and SS, respectively. The propagators are contracted to form baryon 
blocks projected to fixed momentum at the sink for use in the calculation of the correlation 
functions to be described below. The blocks are denned as 

Sg^p,f;*o) =Ee ip-x S^ , (1) 

X 

where is a quark propagator of flavor /, and the indices are combined spin-color indices 
running over i — 1, . . . , N C N S . 1 The choice of the fi and the tensor b^ H > depend on the spin 
and flavor of the baryon, H, under consideration. For our calculations we used the local 
interpolating fields constructed in Ref. [33], restricted to those that contain only upper spin 
components (in the Dirac spinor basis). This choice results in the simplest interpolating fields 
that also have the best overlap with the octet-baryon ground states. Blocks are constructed 
for all lattice momenta |p| 2 < 4 allowing for the study of multibaryon systems with zero or 
nonzero total momentum and with nontrivial spatial wave functions. 

B. Multibaryon Interpolating Operators and Contractions 

In order to define correlation functions for the multihadron systems, interpolating operators 
with well defined quantum numbers at the source and sink are constructed. As we intend 



1 To be specific, for a quark spin component i s — l,...,N a and color component i c — 1,...,N C , the 
combined index % — N c (i s — 1) + i c . 
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to perform calculations away from the SU(3)-flavor symmetry limit at lighter quark masses, 
the quantum numbers of parity 7r, angular momentum J 2 and J z , strangeness s, baryon 
number (atomic number) A, and isospin I 2 and I z are used to define the interpolating oper- 
ators. 2 These interpolating operators are first constructed recursively at the hadronic level 
from the octet-baryon field operators using the appropriate group products (Clebsch-Gordan 
coefficients for isospin and angular momentum) to build an outer product wavefunction 
|space) (g) |ang.mom.) <g> |isospin) <8> (parity) of given strangeness and baryon number. This 
approach is similar to that used in Ref. [33] in the context of excited baryons. The baryons 
within this wavefunction are then replaced by appropriate quark-level wavefunctions, of 
which there are, in principle, multiple choices, and then a quark-level anti-symmetrization 
is performed (as color is included in the quark level wavefunctions). A similar approach has 
been used to investigate the system [34]. 

The quantum numbers defining the systems that we discuss in this paper are shown 
in Table II. States are given a representative hadronic label (first column in Table II) 
indicating one component of their hadronic level wavefunction. In order to determine which 
SU(3) irreducible representations (irreps) are present in the correlation functions, the states 
are acted on by the quadratic and cubic SU(3) Casimir operators, and by V-spin, U-spin 
and isospin raising and lowering operators, the results of which are presented in Table II 
(eigenvalues of the Casimir operator for relevant SU(3) irreps are tabulated in Appendix A). 
Because of the overall antisymmetric nature of allowed quark-level wavefunctions, a number 
of the constructed interpolating operators give rise to correlation functions that contain only 
one SU(3) irrep, while others contain more than one. 

Given the blocks discussed in the previous section and the quark- and hadron-level wave- 
functions introduced previously, the contractions are performed using an algorithm that is 
described in more detail in Ref. [35]. For a given set of quantum numbers, denoted by Q, we 
have a basis of N w f hadron-level and quark-level wavefunctions, and respectively 
for % — 1, . . . , N w{ . Note that N wi depends on Q. In this work the spatial wavefunction at 
the source is restricted to a single point. In addition, the single-baryon interpolating fields 
are restricted to the upper spin components (in the Dirac basis) only. These two restrictions 
drastically reduce both the size of the space of allowed quark-level wavefunctions, and the 
number of terms each wavefunction can have. In all cases, an orthonormal basis of wave- 
functions consistent with the above constraints is obtained. The construction, as well as 
the simplification of the wavefunctions, is done automatically with symbolic manipulation. 
Finally, after the construction of the wavefunctions, independent checks of transformation 
properties of these wavefunctions were performed, confirming that these wavefunctions trans- 
form as expected. As discussed previously, hadron-level wavefunctions and hadronic blocks 
with a given total momentum are used at the sink. These basic building blocks allow for 
the construction of more interpolating fields at the sink with nontrivial spatial hadronic 
wavefunctions. In addition, hadron systems with non-vanishing total momentum can be 
constructed, since the point sources couple to all momenta. 

The contraction algorithm is then straightforward and amounts to selecting the appro- 
priate indices in all possible ways from the hadron blocks building the hadronic-level sink 
wavefunction, dictated by the quark-level wavefunction. 3 For all the systems studied here, 



2 For calculations restricted to the SU(3)-flavor symmetric limit, it would also be advantageous to work 
directly with SU(3) irreducible representations. 

3 We note that the algorithm proposed in Ref. [36] is quite similar to the one we have been using in the 
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Label 


A 


s 


I 


J n 


Local SU(3) irreps 


This work 


N 


1 





1/2 


1/2+ 


8 


8 


A 


1 


-1 





1/2+ 


8 


8 


£ 


1 


-1 


1 


1/2+ 


8 


8 




1 


-2 


1/2 


1/2+ 


8 


8 


d 


2 








1+ 


TO 


TO 


nn 


2 





1 


0+ 


27 


27 


nA 


2 


-1 


1/2 


0+ 


27 


27 


nA 


2 


-1 


1/2 


1+ 


8 A , TO 


— 


nS 


2 


-1 


3/2 


0+ 


27 


27 


nS 


2 


-1 


3/2 


1+ 


10 


10 


nH 


2 


-2 





1+ 


8 A 






2 


-2 


1 


1+ 


8 A , 10, TO 


— 




2 


-2 





0+ 


1, 27 


1, 27 


3 H, 3 He 


3 





1/2 


1/2+ 


35 


35 


iH(l/2+) 


3 


-1 





1/2+ 


35 


— 


A H(3/2+) 


3 


-1 





3/2+ 


TO 


TO 


He^H, nnA 


3 


-1 


1 


1/2+ 


27, 35 


27, 35 


|He 


3 


-1 


1 


3/2+ 


27 


27 


4 He 


4 








0+ 


28 


28 


A He > A H 


4 


-1 


1/2 


0+ 


28 




AlHe 


4 


-2 


1 


0+ 


27, 28 


27, 28 


AE°pnn 


5 


-3 





3/2+ 


TO + ... 


TO 



TABLE II: The baryon number A, strangeness s, total isospin /, total spin and parity J w quantum 
numbers of the states and interpolating operators studied in the current work. For each set of 
quantum numbers, the SU(3) irreps that are possible to construct with local interpolating operators 
are listed. The last column lists the SU(3) irrep(s) of the interpolating operators used in this work, 
and the dashes indicate that the state is inferred from other states using SU(3) symmetry. 

the total contraction time was an order of magnitude less than the rest of the calculation. 
In addition, the biggest contraction burden was due to the large number of terms contribut- 
ing to the wavefunctions with a nontrivial spatial part at the sink (moving hadrons at the 
sink). As an example of the speed of our contraction code, a 4 He correlation function can 
be computed in ~ 0.8 s per time-slice on a single core of a dual core AMD Opteron 285 
processor. 



production of the results presented here. 
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TABLE III: The pion energy (l.u.) as a function of momentum (l.u.), |P| = (xj l n l) calculated on 
each ensemble of gauge-field configurations. The infinite-volume pion mass, determined by fitting 
the expression in eq. (2), is provided in the last row. The first uncertainty is statistical and the 
second is the fitting systematic. 



ensemble 


|n|=0 [n| 2 = 1 |n| 2 = 2 |n| 2 = 3 |n| 2 = 4 |n| 2 = 5 


24 3 x 48 
32 3 x 48 
48 3 x 64 


0.59389(18) (18) 0.64652(16) (19) 0.69482(17)(29) 0.73971(20)(36) 0.77800(30) (72) 0.81946(36) (78) 
0.59445(15) (17) 0.62474(15) (18) 0.65326(16)(20) 0.68099(18)(25) 0.70672(19)(28) 0.73194(22) (31) 
0.59403(16) (14) 0.60768(16)(15) 0.62101(18)(17) 0.63403(19)(20) 0.64667(21)(24) 0.65915(24)(28) 


L = oc 


0.59426(12) (11) 



III. THE PION AND BARYON DISPERSION RELATIONS 



In the limit of SU(3)-flavor symmetry, all members of the lightest baryon octet have the 
same mass, and as such, we compute correlation functions associated with only one of 
the octet baryons. Similarly all octet pseudoscalar mesons are degenerate, and we refer to 
them as the pion. Linear combinations of single-hadron correlation functions generated from 
smeared quark sources and either smeared or point sinks are formed for hadrons with a given 
lattice momentum. The lowest energy eigenvalue can be determined from these correlation 
functions, the results of which are presented in Table III (pion) and Table IV (baryon), and 
the baryon effective mass plots (EMPs) are shown in fig. 1. 



■•••'-*T**tj| 



4 8 12 16 20 4 8 12 16 20 4 8 12 16 20 

t/b t/b t/b 





FIG. 1: The EMPs associated with linear combinations of baryon correlation functions computed 
with the 24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles, with momentum 
|P| = 0. The inner (darker) shaded region corresponds to the statistical uncertainty of the extracted 
energy, while the outer (lighter) shaded region corresponds to the statistical and fitting systematic 
uncertainties combined in quadrature. The time-extent of each band corresponds to the choice of 
the fitting interval for each correlation function. 

For hadrons at rest, the masses of the pion and baryon in finite volume, m^^m,!), are 
extrapolated to infinite volume using 

MP(m n L) = MjT + 4 y) e —— + ... , (2) 

where only the first terms in the finite-volume (FV) expansion are required due to the 
large pion mass [37]. The extrapolations to infinite volume are shown as the solid regions 
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TABLE IV: The ground-state octet-baryon energy (l.u.) as a function of momentum (l.u.), 
|P| = (77) l n l' calculated on each ensemble of gauge-field configurations. The infinite-volume 
baryon mass, determined by fitting the expression in eq. (2), is provided in the last row. The first 
uncertainty is statistical and the second is the fitting systematic. 



ensemble 


|n| = 


|n| 2 = l |n| 2 = 2 |n| 2 = 3 |n| 2 = 4 


|n| 2 = 5 


24 3 x 48 
32 3 x 48 
48 3 x 64 


1.20317(58) (84) 
1.20396(47)(69) 
1.2032(07)(11) 


1.2282(9)(16) 1.2537(9)(23) 1.2785(11)(31) 1.3023(11) (25) 
1.21821(61)(64) 1.23263(65)(70) 1.24685(69)(79) 1.26077(74) (94) 
1.2096(11) (22) 1.2162(11) (21) 1.2227(12) (22) 1.2290(12) (21) 


1.3254(12) (29) 
1.2746(08) (11) 
1.2354(13) (21) 


L = oo 


1.20293(35) (20) 



in fig. 2, and the extrapolated values of the pion and octet-baryon mass are presented 
in Table III and Table IV, respectively. As expected for calculations with large values 
of m^L, the single-hadron FV effects are very small. The extrapolated pion and octet- 



0.5950 
0.5945 
£ 0.5940 
0.5935 



0.5930 
10 




10" 



10" 



10" 



(m n L) 



i/2 



Exp[- m K L] 



ur 




{m„L) 



FIG. 2: The volume dependence of the pion mass (left panel) and the baryon mass (right panel) 
extracted from the zero-momentum correlation functions. The shaded regions are extrapolations 
of the form given in eq. (2). 

baryon masses, using the measured lattice spacing, are m n = 805.9(0.6) (0.4) (8.9) MeV and 
tub = 1.634(0) (0) (18) GeV, where the first uncertainty is statistical, the second is the fitting 
systematic, and the third is due to the uncertainty in the lattice spacing. 

In order to have confidence in the extraction of multibaryon binding energies and to be 
able to quantify one of the systematic uncertainties in these determinations, it is important 
to determine the single-hadron dispersion relation. The energies of the pion and baryon are 
shown in fig. 3 as a function of J2 sin 2 ( ^nA, where the triplet of integers n = (72i,n 2 ,n 3 ) 

3 

is related to the lattice momentum via |P| 2 = f^M |n| 2 . In these LQCD calculations, the 
energy of the hadron can related to its lattice momentum through a dispersion relation of 
the form 

(bE H f = (bM H f + -L £ sin 2 (^n^j , (3) 
where the anisotropy parameter, (or equivalently the speed of light c = !/£#), is expected 
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j j 

FIG. 3: The squared energy (in (l.u.) 2 ) of the single pion and baryon as a function of 
The points are the results of the LQCD calculations with the inner (outer) 

j 

uncertainties being the statistical uncertainties (statistical and systematic uncertainties combined 
in quadrature). The red curves correspond to the best linear fits. 

to be unity in calculations performed with isotropic lattices. 4 Fitting to the energy of the 
pion and baryon, given in Table III and Table IV, respectively, yields £ n = 1.0055(57) (26) 
and = 1.019(10)(03). Therefore, the dispersion relations provide only a small uncertainty 
in the extraction of multihadron energies. 

IV. TWO-BODY SYSTEMS 

In general, the two-body states can be classified by isospin, strangeness, parity and angular 
momentum. In the limit of SU(3)-flavor symmetry, the energy eigenstates can also be 
classified by SU(3) quantum numbers. The lowest-lying baryons transform as 8 under SU(3), 
and, therefore, the two-body states have degeneracies determined by the dimensionality of 
the irreps in the product 

8©8 = 27©10©T0©8 5 ffi8 A ©l. (4) 

As the wavefunction of such systems is totally antisymmetric, the s-wave channels trans- 
form under SU(3) as 27©8g©l, while the 3 Si- 3 Di coupled channels transform as 10©10©8 A . 
The source structures we have employed, in which the quark-level operators reside at one 
point in the spatial volume, have vanishing overlap with the 85 irrep, and as a result, we 
are unable to determine the energy of this two-body irrep. Correlation functions are not 
constructed directly in terms of their SU(3) transformation properties, but the contributing 
SU(3) irreps can be deduced from their structure: 10 from the deuteron, 27 from the di- 
nucleon, 1 ©27 from the H-dibaryon (the 85 is absent), 10 from nS~ in the 3 S\- 3 Di coupled 
channels, and 8 A from 1 = iVS in the 3 S±- 3 Di coupled channels. EMPs extracted from the 
two-body correlation functions for systems at rest calculated with the 48 3 x 64 ensemble 



4 As the lattice hadronic dispersion relations are a priori unknown, they must be calculated. The form 
given in eq. (3) is expected to capture the leading momentum dependence. 
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FIG. 4: EMPs associated with |P| = two-baryon correlation functions computed with the 48 3 x 
64 ensemble. The inner (darker) shaded region corresponds to the statistical uncertainty of the 
extracted energy, while the outer (lighter) shaded region corresponds to the statistical and fitting 
systematic uncertainties combined in quadrature. The time-extent of each band corresponds 
to the choice of the fitting interval for each correlation function. From left to right, the top 
row corresponds to the 1, 8,4, 10 SU(3) irreps (corresponding to the H-dibaryon, 1 = NE in the 
3 5i— 3 Di coupled channels and n£~ in the 3 S\ — 3 D\ coupled channels, respectively), and the bottom 
row corresponds to the 10 and 27 (corresponding to the deuteron and di-neutron, respectively). 



are shown in fig. 4. The energies of states that are negatively shifted relative to two free 
baryons are presented in Table V, Table VI and Table VII, respectively, and displayed in 
fig. 5. 

The energies of the states that are presented in this work, along with their statistical 
uncertainties, are determined from a single-parameter correlated ^-minimization procedure 
performed over a specific time interval of EMPs and from exponential fits to the correlation 
functions directly, with covariance matrices determined with either Jackknife or Bootstrap. 
The systematic uncertainty that is assigned to these energies is determined by varying the 
fit interval over a range of values consistent with the identified plateau region. 

A number of scattering states with positive energy-shifts relative to two free baryons have 
also been identified using different correlation functions, but their uncertainties are large 
enough to preclude clean extraction of scattering phase-shifts using Luschers method [38, 39], 
and we defer analysis of these states to a later time when adequate statistics have been 
accumulated. 

In sufficiently large volumes, the binding momentum associated with a two-body bound 
state at rest in the lattice volume will scale as 

k(L) = k + -^e" Koi + ... , (5) 

where kq is the infinite-volume binding momentum, k,q = VMbB, where B is the binding 
energy and is the residue of the bound-state pole [38-40]. Analogous FV scaling formulas 
for systems moving in the lattice volume are known [26], but at this order in the expansion 
differ from the relation in eq. (5) only by the coefficient of the second term. In the 32 3 x 48 
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TABLE V: Two-body binding energies (MeV) calculated with the 24 3 x 48 ensemble. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. 



SU(3) irrep 


|n| = 


|n| = 1 


|n| = 2 


1 


77.7(1.8)(3.2)(0.8) 


67.2(2.5)(2.5)(0.8) 


85.0(3.1)(4.0)(0.9) 


8a 


40.1(1.7)(2.9)(0.4) 


26.5(1.8)(3.6)(0.3) 


46.7(2.0)(3.2)(0.5) 


10 


11.4(1.8)(4.0)(0.1) 


6.3(1.9)(4.4)(0.1) 


15.3(2.2)(4.5)(0.1) 


TO 


25.4(2.6)(4.7)(0.3) 


16.0(2.7)(5.9)(0.2) 


40.7(3.6)(7.4)(0.5) 


27 


17.8(1. 7)(2.8)(0.2) 


6.9(1.8)(3.8)(0.1) 


28.5(2.3)(3.8)(0.3) 



TABLE VI: Two-body binding energies (MeV) calculated with the 32 3 x 48 ensemble. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. 



SU(3) irrep 


|n| = 


|n| = 1 


|n| = 2 


1 


76.0(2.3)(2.8)(0.8) 


70.3(2.3)(3.1)(0.7) 


79.6(2.6)(3.9)(0.9) 


8a 


38.5(2.3)(4.4)(0.4) 


34.0(2.6)(3.4)(0.4) 


45.2(3.0)(3.1)(0.5) 


10 


10.5(2.5)(4.1)(0.1) 


1.1(2.4)(4.2)(0.0) 


12.9(2.6)(4.5)(0.1) 


TO 


22.5(2.3)(2.6)(0.2) 


19.2(2.3)(3.7)(0.2) 


31.6(2.7)(3.2)(0.3) 


27 


15.1(2.0)(2.0)(0.2) 


12.3(1.9)(3.6)(0.1) 


24.9(2.2)(3.1)(0.3) 



and 48 3 x 64 lattice volumes, the energies of the two-body bound states do not exhibit 
statistically significant volume dependence. Consequently, using eq. (5) to determine the 
infinite- volume binding energies does not provide a refinement over simply taking the binding 
energies determined in the 48 3 x 64 ensemble, and the latter is used as the best estimate 
of the infinite-volume binding energies, the results of which are shown in Table VII. The 
expected differences between the infinite- volume bindings and those in the 48 3 x 64 ensemble 
can be estimated from the values of k L given in Table VII. With the exception of the state 

TABLE VII: Two-body binding energies (MeV) calculated with the 48 3 x 64 ensemble. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. The second to last column corresponds to an average of the |n| = 0,1,2 calculations, 
which is taken to be the infinite volume value. The last column gives the value of kq times the 
spatial lattice size for L = 48. 



SU(3) irrep 


|n| = 


|n| = 1 


|n| = 2 


L = oo 


KqL 


1 


73.7(3.3)(5.1)(0.8) 


73.7(4.4)(7.6)(0.8) 


75.4(3.3)(3.3)(0.8) 


74.6(3.3)(3.3)(0.8) 


12.3 


8a 


38.7(2.9)(2.9)(0.4) 


34.6(2.8)(3.1)(0.4) 


39.7(3.0)(2.7)(0.4) 


37.7(3.0)(2.7)(0.4) 


8.8 


10 


6.6(3.4)(4.1)(0.0) 


2.8(3.1)(4.1)(0.0) 


7.0(3.4)(3.7)(0.0) 


5.5(3.4)(3.7)(0.0) 


3.3 


TO 


19.7(3.1)(4.1)(0.2) 


17.8(3.6)(3.1)(0.2) 


23.1(3.9)(5.5)(0.2) 


19.5(3.6)(3.1)(0.2) 


6.3 


27 


13.1(2.8)(4.3)(0.2) 


14.9(2.7)(2.7)(0.2) 


19.3(2.9)(3.3)(0.2) 


15.9(2.7)(2.7)(0.2) 


5.7 
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FIG. 5: Binding energies in the A = 2 systems relative to two non-interacting baryons (B = 
—AE). The points and associated uncertainties are the results of the LQCD calculations given 
in Table V, Table VI and Table VII. The dark (statistical uncertainty) and light (statistical 
and systematic uncertainties combined in quadrature) horizontal bands denote the average of the 
bindings calculated on the 48 3 x 64 ensemble, which are taken as the infinite-volume estimate. 
Where only bands are shown, SU(3) symmetry has been used to determine the spectrum. 

in the 10 irrep, the states are small enough compared to the lattice volume to make the 
finite-volume effects negligible. 

There are a few important results that should be highlighted. The deuteron is 
found to be substantially more deeply bound in the present calculations, B^ f ~ 3 = 
19.5(3.6)(3.1)(0.2) MeV, than in the quenched calculations [8] in which a binding energy 
of B n d f = 9.1(1. 1)(0.5) MeV at a similar pion mass is found. The H-dibaryon is found 
to be deeply bound with B H = 74.6(3.3) (3.3) (0.8) MeV, approximately twice as bound as 
the result found by HALQCD [16] at a similar quark mass. In recent work we reported 
that the interaction in the 3 Si — 3 Di channel was extremely repulsive at a pion mass 
of Tti-ji ~ 390 MeV [23], consistent with the phase-shift analysis of experimental data at the 
physical pion mass. At the SU(3) symmetric point, we find that this state has moved close to 
threshold and is even consistent with being bound, indicating that there is significant light 
quark mass dependence in this channel at the heavier quark masses (beyond the regime of 
applicability of the relevant EFT). 

As the calculations have been performed at the SU(3) symmetric point, the states dis- 
cussed above provide a nearly complete set of two-baryon ground states, with only the 8s 
irrep being absent. Furthermore, since the determinations of the various energy levels in 
the two-body sector are correlated, their differences can be determined more precisely than 
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TABLE VIII: Two-body energy splittings, Ej 1 — Ej 2 (MeV) between different multiplets calculated 
with the 32 3 x 48 ensemble. The column refers to representation I\ and the row to representation 
h- The first uncertainty is statistical, the second is the fitting systematic and the third is due to 
the lattice spacing. 



h\h 


S A 


10 


TO 


27 


1 


34.3(0.7)(1.2)(0.4) 


65.9(0.4)(0.9)(0.7) 


49.3(1.4)(1.7)(0.5) 


55.4(1.2)(1.8)(0.6) 


8 A 




31.0(0.8)(1.6)(0.3) 


14.2(1. 1)(2.0)(0.2) 


20.5(1.2)(2.3)(0.2) 


10 






-17.7(1.4)(2.2)(0.2) 


-11.0(1.4)(2.4)(0.1) 


TO 








5.8(1.0)(1.0)(0.1) 



their individual values. In Table VIII we present the splittings between the various irreps. 
The energy difference between the 10 and 27 corresponds to the deuteron-di-nucleon mass 
difference. This splitting is found to be small, and consistent with zero within the uncer- 
tainties of the calculation. Theoretically, it has been established from SU(2) that these 
states become degenerate in the large- N c limit of QCD [41], with a fractional splitting (and 
violation of Wigner's SU(4) symmetry) that scales as 1/iV^. Extending the argument to 
the strange sectors shows that the other splittings are only l/N c suppressed, and not 1/N% 
suppressed [41]. Such scalings are consistent with what we have found, but verification of 
the scaling will require significantly higher statistics in the calculations. 

V. THREE-BODY SYSTEMS 

The correlation functions for the three-body systems are generated using the procedure 
described previously. As is the case for two-body systems, the states in the spectrum for 
each system can be classified by their SU(3) quantum numbers in the limit of SU(3)-flavor 
symmetry. The three-body states can be assigned to the SU(3) irreps in 8 ® 8 <g> 8, which 
can be straightforwardly constructed as 

8©8©8 = 64©2 35©235©6 27©4 10©4T0©88©21 . (6) 

However, the local sources constructed from only the upper components of the quark fields 
produce correlation functions containing a subset of these irreps, 

8©8©8^35©35©2 27©10©T0©28©1 , (7) 

and further decomposition into states with J w = \ + and J w = | + gives 

( 8 © 8 © 8 ) JW=1/2+ ->■ 35 © 35 © 27 © 8 

( 8©8©8 ) JW=3/2+ ->■ 27© 10© T0©8© 1 . (8) 

It is clear from the SU(3) irreps contributing to the three-body systems that, with our 
source structure, a given correlation function contains contributions from multiple SU(3) 
irreps. With a relatively small number of states identified with the present set of correlation 
functions, the SU(3) classification of states is difficult to establish from the spectra alone. 
More generally, it is expected that the spectrum of states in any given correlation function 
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becomes increasingly complicated with increasing numbers of baryons even when constrained 
by SU(3)-flavor symmetry. As the focus of this work is systems containing only a small 
number of strange quarks, we have chosen to use the same notation as in hypernuclear 
spectroscopy. States in 3 He (same as 3 H by isospin symmetry), ^He (same as and nnA 
by isospin symmetry), the isosinglet ^H, and the isotriplet |He have been identified in the 
three-body sector. 

Correlation functions calculated with LQCD will not only contain contributions from the 
ground state and excited states of the bound nuclei but also continuum states that consist 
of all possible sub-clusterings of the baryons. For instance, the correlation functions used 
to extract the 3 He nuclear states will also contain contributions from the deuteron-proton 
and di-proton-neutron in addition to the proton-proton-neutron continuum states. With 
sufficient precision in the calculation, one will be able to use these levels to extract, for 
instance, the deuteron-proton scattering phase-shift [25] . Given that the two-body sector is 
well-established, the spectrum of such continuum states can be approximately constructed. 
Clearly, states of the 3 He nucleus can only be cleanly identified when they are not close in 
energy to the expected location of non-interacting continuum states. The generalization of 
this discussion applies to other systems comprised of three or more baryons. In Appendix 
B, an example of the expected FV scattering-state spectrum is constructed for each of the 
volumes used in this analysis, demonstrating the extent of this problem in large volumes. 



1 jn 
2' J 



3 H and 3 He 



In nature, the I = |, J T = | + ground state of the 3 He nucleus is the only bound state of two 
protons and a neutron, and it is known to be dominantly composed of two protons in a 
state coupled to an s-wave neutron. Four 3 He correlation functions, resulting from different 
source structures defined by s = 0, / = | and J n = ^ + quantum numbers transforming 
as a 35 of SU(3), have been constructed. 5 EMPs obtained from correlation functions in 
each of the three ensembles, from which the energy of the lowest-lying 3 He states have been 
determined, are shown in fig. 6. 



0.10 
0.05 

W 

< o 

-0.05 
-0.10 



8 12 
t/b 



0.10 
0.05 

W 

< o 
£> 

-0.05 
-0.10 



8 12 
t/b 



0.10 
0.05 

W 

< o 
X! 

-0.05 
-0.10 



"-7T1 



8 12 
t/b 



FIG. 6: EMPs associated with J n = 1 + 3 He ( 3 H) |P| = correlation functions computed with the 
24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles. The inner (darker) shaded region 
corresponds to the statistical uncertainty of the extracted energy, while the outer (lighter) shaded 
region corresponds to the statistical and fitting systematic uncertainties combined in quadrature. 



The only possible SU(3) irrep with these quantum numbers is the 35. 
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TABLE IX: The calculated J 71 " = \ binding energy of 3 He ( 3 H) in the 24 3 x 48 ensemble, "g.s." 
denotes the ground state. The first uncertainty is statistical, the second is the fitting systematic 
and the third is due to the lattice spacing. 



3 He 


24 3 x 48 




|n| 2 = 


|n| 2 = l 


|n| 2 = 2 


g.s. (MeV) 


65.4(5.1)(4.4)(0.7) 


42.8(3.8)(8.9)(0.4) 


46.3(5.3)(6.7)(0.5) 



TABLE X: The calculated J w = \ binding energy of 3 He ( 3 H) in the 32 3 x 48 ensemble, "g.s." 
denotes the ground state. The first uncertainty is statistical, the second is the fitting systematic 
and the third is due to the lattice spacing. 



3 He 


32 3 x 48 




|n| 2 = |n| 2 = 1 


|n| 2 = 2 


g.s. (MeV) 


63.2(3.9)(7.0)(0.7) 52.9(5.7)(9.9)(0.6) 


55.7(6.4)(10.1)(0.6) 



The 3 He bound-state energies on the ensembles are given in Tables IX, X and XI and are 
shown in fig. 7 along with the thresholds for non-interacting continuum states. 6 The exact 
form of infinite- volume extrapolation of three- and higher-body bound-state energies is as yet 
unknown, though expected to be exponential (see Refs. [42-44] for related discussions). For 
the current study, we simply average the results obtained from the system at rest and from 
the boosted systems on the 48 3 x 64 ensemble to provide an estimate of the infinite-volume 
binding energy of 

#(°o)(3 He ) = 53.9(7.i)(8.0)(0.6) MeV . (9) 

The energy of this state is significantly lower than any of the expected continuum states, 
based upon where they would lie in the spectrum in the absence of interactions. Therefore, 
we conclude that this is the ground state of 3 He. 

While it is tempting to compare these results with the experimental spectrum of 3 He, 
one should refrain at present, since these calculations are performed in the SU(3) limit of 



TABLE XI: The calculated J 71 " = \ binding energy of 3 He ( 3 H) in the 48 3 x 64 ensemble, "g.s." 
denotes the ground state. The first uncertainty is statistical, the second is the fitting systematic 
and the third is due to the lattice spacing. 



3 He 


48 3 x 64 






n 


2 = 


n 


2 = 1 


n 


2 = 2 


g.s. (MeV) 


61.9(8 


.9)(10.9)(0.7) 53.0(7.1)(8.0)(0.6) 50.0(6.1)(9.2)(0.6) 



Finite-volume effects will lead to small shifts in these thresholds. 
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FIG. 7: The bound-state energy levels in the J 71 " = ^ 3 He ( 3 H) sector. The points and their 
associated uncertainties correspond to the energies of the states extracted from the correlation 
functions with the quantum numbers of the ground state of 3 He. The locations of the scattering 
thresholds associated with non-interacting deuteron-proton, di-proton-neutron and proton-proton- 
neutron continuum states, determined from the single-hadron spectrum and the two-body binding 
energies given in Table VII are shown. 

QCD and without electromagnetism. The ground-state binding energy will receive a shift 
due to the electromagnetic interaction between the two protons. On the other-hand, the 
exact isospin symmetry directly relates this spectrum to that of the triton. In nature the 
triton binding energy per nucleon is B/A ~ 2.83 MeV, while at the SU(3) symmetric point 
we find that B/A ~ 24 MeV, more than an order of magnitude larger. 

The 3 He ground-state energy that we have calculated in this rif = 3 calculation is substan- 
tially different from that obtained with quenched calculations at a comparable pion mass [9] , 
which find an infinite- volume extrapolated value of 2?^°2 ( 3 He) = 18.2(3.5) (2.9) MeV. A 
likely explanation for the difference is quenching artifacts, which are unlikely to cancel be- 
tween the bound system and the threshold states. The difference in the total energy (not the 
binding energy) of the 3 He ground state between the two calculations is of 0(1%), smaller 
than the differences observed between single-hadron masses in quenched and unquenched 
calculations [45]. Additionally, the contributions from continuum states that must be present 
in both calculations at some level (see Appendix B) may pollute the extraction of the 3 He 
ground state, particularly in large volumes. 
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TABLE XII: The calculated J 77 = § binding energies in ^H. "g.s." denotes the ground state. The 
energies in the J n = ^ + channel are the same as those of 3 He by SU(3) symmetry, see Tables IX, 
X and XI. The first uncertainty is statistical, the second is the fitting systematic and the third is 
due to the lattice spacing. 



in 


24 3 x 48 


32 3 x 48 


48 3 x 64 


J w = l + g.s. (MeV) 


90.8(4.5)(6.5)(1.0) 


89.6(4.6)(8.9)(1.0) 


82(8)(12)(1) 



B. 



0, J 77 



i + 

2 



and J n 



3 H 



The Hypertriton 



The hypertriton, ^H, with the quantum numbers of npA and / = is the simplest hyper- 
nucleus produced in the laboratory, having a total binding energy of B ~ 2.35 MeV. With 
a A-separation energy of just B k ~ 0.13 MeV, it is consistent with a A weakly bound to a 
deuteron. The ground state has J n = | + and has been identified as a member of the 35 of 
flavor SU(3) [46]. It continues to be the focus of experimental efforts, for instance, in heavy- 
ion collisions at RHIC [47] and the HypHI project at GSI, where in the latter it is being 
used as a "phase-zero" calibration nucleus for the production and detection systems [48] . We 
have calculated correlation functions in both the J n = | + and J 71 " = | + channels and have 
identified the lowest-lying state in each. Two of the correlation functions associated with 
the J n = \ + channel are pure 35 and are in the same irrep as 3 He, and hence the energy 



of the identified states are the same. Further, the J 71 " 
the J 77 



| + channel is pure 10. EMPs in 

o _|_ 

| channel from these correlation functions are shown in fig. 8, from which the 
energies of the lowest lying states have been determined, and are given in Table XII. The 
EMPs in the J n = | + channel are not shown, as they are identical to those of 3 He, shown 
in fig. 6. The extracted spectra of bound states are shown in fig. 9. Taking the results 
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FIG. 8: The EMPs associated with a J 71 " = | + hypertriton (^H) correlation function computed with 
the 24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles, with momentum |P| = 0. The 
inner (darker) shaded region corresponds to the statistical uncertainty of the extracted energy, while 
the outer (lighter) shaded region corresponds to the statistical and fitting systematic uncertainties 
combined in quadrature. 

obtained in the 48 3 x 64 ensemble to be the best estimate of the infinite-volume binding 
energies gives, 



J B (oo) ( 3 v H(l/2 + )) = 53.9(7.1)(8.0)(0.6) MeV 
J B (oo) (lH(3/2+)) = 82(8)(12)(1) MeV , 



(10) 
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FIG. 9: The bound-state energy levels in the J n = \ (upper panel) and J 7r = | (lower panel) 
hypertriton (^H) sector. The points and their associated uncertainties correspond to the energies of 
the states extracted from the correlation functions with the quantum numbers of the ground state 
of ,r = \ and ,r = | |H. The locations of the energy-levels associated with non-interacting 
continuum states, determined from the two-body binding energies given in Table VII, are shown. 
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where we have used the 3 He result for the J 77 = | binding energy, which includes calcula- 
tions of boosted systems. 

The observed states are significantly below the scattering thresholds and are consistent 
with a bound nucleus at these values of the quark masses in the absence of electro- 
magnetism. Interestingly, the lowest energy state is in the J 71 " = | + spin channel. As the 
measurements of the two spin states are correlated, the spin splitting can be extracted with 
high precision, resulting in 

fi (oo) (lH(3/2+)) - fi (oo) (iH(l/2+)) = 26.2(2.3)(5.5)(0.3) MeV . (11) 



1, J n 



jjie, |H and nnA 



The isotriplet of states 7 , ^He, and nnA, are degenerate in the absence of electromag- 
netism and in the limit of exact isospin symmetry, and can have J n = \ and J 77 = § . The 
J n = | + is expected to be the lowest-lying state, with a significant component of the wave- 

function having the two nucleons in the So channel coupled to the A. The J w = | state 
cannot have such a NNA configuration in its wavefunction by the Pauli principle without 
placing the baryons in orbital excitations but will have configurations of the form of two 
nucleons in the 3 Si- 3 Di channel coupled to a S + . In the SU(3) limit, this can be nearby 
in energy, but when SU(3) breaking is included, the energy for the J 77 = | will increase, 
largely dictated by the E-A mass splitting, and become less phenomenologically interesting. 
Consequently, we will focus first on the J 71 " = | + channel. The EMPs from one of the eight 
correlation functions of these quantum numbers are shown in fig. 10, from which the energies 
of the lowest-lying states have been determined. The extracted spectrum of bound states 




^He (^H and nnA) correlation function computed 



FIG. 10: The EMPs associated with one J 71 " = 
with the 24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles, with momentum 
|P| =0. The inner (darker) shaded region corresponds to the statistical uncertainty of the extracted 
energy, while the outer (lighter) shaded region corresponds to the statistical and fitting systematic 
uncertainties combined in quadrature. 



is given in Table XIII, and shown in fig. 11. Taking the result obtained on the 48 3 x 64 
ensemble as the estimate of the infinite-volume binding energy, we find 



5 (oo) (iHe(l/2+)) = 69(5)(12)(0) MeV 



(12) 



7 We refer to the npA state with the np coupled to / 
which the np couple to / = 0. 



1 as to differentiate it from the state in 
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TABLE XIII: The calculated binding energies in ^He (^H and nnA). "g.s." denotes the ground 
state. The first uncertainty is statistical, the second is the fitting systematic and the third is due 
to the lattice spacing. 



iHe 


24 3 x 48 


32 3 x 48 


48 3 x 64 


g.s. (MeV) 


77.6(3.6)(7.5)(0.8) 


74.1(3.9)(7.3)(0.8) 


69(5)(12)(0) 
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FIG. 11: The bound-state energy levels in the J n = \ ^He (^H and nnA) sector. The points and 
their associated uncertainties correspond to the energies of the states extracted from the correlation 
functions with the quantum numbers of the ground state of }\He. The locations of the energy-levels 
associated with non-interacting di-proton-A, AN-N and A-N-N continuum states, determined from 
the two-body binding energies given in Table VII, are shown. 

The ground state is significantly more deeply bound than any of the continuum states, and 
we identify this as the ground state of the ^He nucleus (and hence also bound and nnA 
due to isospin symmetry). The correlation function from which this ground-state energy 
was extracted is a superposition of 35 and 27 SU(3) irreps. Another element of the 27 
irrep is in the s = —3 sector, with I = 1, J n — | and with the baryon structure of NEA. 
One of the correlation functions associated with this state is pure 27, and the energy of the 
lowest-lying state in this correlation function is found to be the same as that in the ^He 
correlation function within the uncertainties of the calculations, suggesting that the 27 state 
is lower in energy than or nearly degenerate with the 35. 

Experimentally, there is no evidence for a bound ^He nucleus as the A-nucleon interactions 
are not sufficient to overcome the Coulomb repulsion between the protons. Further, the small 
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TABLE XIV: The calculated binding energies in J n = | fjHe. "g.s." denotes the ground state. 
The first uncertainty is statistical, the second is the fitting systematic and the third is due to the 
lattice spacing. 



|He 


24 3 x 48 


32 3 x 48 


48 3 x 64 


g.s. (MeV) 


64.3(4.5) (7.9) (0.7) 


58.2(5.2) (7.7) (0.6) 


55(6)(10)(1) 



binding of the hypertriton, with a significant deuteron-A component, strongly suggests that 
the corresponding 1=1 state will be unbound, and it is likely, but yet to be verified, that 
the nnA electrically neutral nucleus is also unbound. However, our calculations provide 
compelling evidence for a bound state in this channel in the limit of SU(3)-flavor symmetry, 
and we expect that the bound state persists over a range of light-quark masses. 



D. 



i, ,r 



|He 



As discussed above, for the 1 = 1, s = — 1, J 7r = | + , three body state, an NNA component 
is forbidden (for all baryons in a relative s-wave) and one important contribution to the 
ground-state wavefunction is pnS, where the nucleons couple to / = 0, J = 1, as in ^H. As 
yet, the only observed S hypernucleus is |He (ppnT, ) [49, 50], but at the SU(3) point it is 
possible that this three-body system binds. The sources used to generate this correlation 
function transform as 27 under SU(3) 8 , and result in EMPs that exhibit clear plateaus. 
The ground-state energies extracted from the three ensembles are given in Table XIV, and 
the associated EMPs are shown in fig. 12. The ground-state energy and the anticipated 
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FIG. 12: The EMPs associated with one J n = | f^He correlation function computed with the 
24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles, with momentum |P| = 0. The 
inner (darker) shaded region corresponds to the statistical uncertainty of the extracted energy, while 
the outer (lighter) shaded region corresponds to the statistical and fitting systematic uncertainties 
combined in quadrature. 

continuum thresholds based upon the non-interacting two-body energies are shown in fig. 13. 



This 27 irrep is different from that in the J x = ^ + channel. In principle the ground state of the system 
could reside in the 64 irrep, but this is not accessible with our present operator structure. 



21 



50 







< 



-50 




























100 



He<| , 



2 + n + p 



d + 2 



• 


L= 


=24 


|p|=0 


• 


L= 


=32 


|p|=0 


• 


L= 


=48 


|p|=0 



EN(3sl) + N 



FIG. 13: The bound-state energy levels in the J n = | |^He sector. The points and their associated 
uncertainties correspond to the energies of the states extracted from the correlation functions with 
the quantum numbers of the ground state of f^He. The locations of the energy-levels associated 
with non- interacting continuum states, determined from the two-body binding energies given in 
Table VII, are shown. 



VI. FOUR-BODY SYSTEMS 



There are a large number of four-body systems and states that could be explored theoretically 
with LQCD at the SU(3) symmetric point, dictated by the product of four 8's, 

8 © 8 © 8 © 8 = 8 1 © 32 8 © 20 10 © 20 10 © 33 27 © 2 28 © 2 28 © 15 35 © 15 35 

© 12 64 ©3 81 ©3 81 ©125 , (13) 

giving a total of 166 lowest-lying states (one per distinct irrep) with distinguishable quantum 
numbers. The local sources that have been used in this work to generate correlation functions 
project onto a subset of the irreps, 

( 8©8©8©8 ) Jn=0+ -»■ 1©27©28 

( 8(8)8® 8(8)8 ) JW=1+ -»■ 8 ©10 ©10 ©35 

( 8 © 8 © 8 © 8 ) Jn=2+ ->■ 8 © 27 , (14) 

which greatly reduces the complexity of individual correlation functions. In order to restrict 
ourselves to systems that are currently of phenomenological importance, we explore systems 
containing up to two strange quarks only, the isosinglet 4 He, the isodoublet A H and \H.e, 
the isosinglet A ^H and the isotriplet AA.He, aaH, and nnAA. 
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TABLE XV: The calculated binding energies in 4 He. "g.s." denotes the ground state. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. 



4 He 


24 3 x 48 


32 3 x 48 


48 3 x 64 


g.s. (MeV) 


115(11)(20)(1) 


107(15)(20)(1) 


107(12)(21)(1) 



a. 1 = 0, r 



H 



In nature, the 4 He nucleus is anomalously deeply bound when compared to nuclei nearby 
in the periodic table, due to its closed shell structure, with a total binding energy of B ~ 
28 MeV, or a binding energy per nucleon of B/A ~ 7 MeV. We anticipate that at the 
SU(3) symmetric point, the binding energy of 4 He will be even deeper given the bindings 
of the deuteron and di-neutron found in the two-body sector. Two of the 4 He correlation 
functions, resulting from different source structures defined by s = 0, I = and J 71 = + 
quantum numbers, transform as an element of the 28 irrep of SU(3) 9 . EMPs of one of these 
correlation functions are shown in fig. 14, from which the energies of the lowest lying states 
have been determined. The extracted spectrum of bound states, only calculated for the 
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FIG. 14: EMPs associated with a |P| = J 71 " = + 4 He correlation function computed with the 
24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles. The inner (darker) shaded region 
corresponds to the statistical uncertainty of the extracted energy, while the outer (lighter) shaded 
region corresponds to the statistical and fitting systematic uncertainties combined in quadrature. 



system at rest in the lattice volume, is given in Table XV and shown in fig. 15. Also shown 
in fig. 15 are the thresholds of non-interacting continuum states, based upon the two-body 
and three-body bound-state spectra. Using the result obtained on the 48 3 x 64 ensemble 
as an estimate of the binding energy in infinite volume, we find 



#M( 4 He) = 107(12)(21)(1) MeV 



(15) 



While this state is somewhat more deeply bound than any continuum state, the precision 
of the calculation is not sufficient to unambiguously distinguish the state from the n + 3 He 
continuum. In order to eliminate this ambiguity in state identification, further calculations 
are required, and additional source structure should be used to increase the size of the basis 
of correlation functions. 



9 The 28 is the only allowed I = 0, s = 0, A=4 irrep. 
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FIG. 15: The bound-state energy levels in the J w = + 4 He sector. The points and their associated 
uncertainties correspond to the energies of the states extracted from the correlation functions with 
the quantum numbers of the ground state of 4 He. The locations of the energy-levels associated with 
non-interacting N- 3 He, d-d, di-nucleon-di-nucleon, di-nucleon-N-N, d-N-N and N-N-N-N continuum 
states, determined from the two-body binding energies given in Table VII and the three-body 
energies given in eq. (9), are shown. 



The 4 He ground-state energy that we have calculated in this rif — 3 calculation is substan- 
tially different from that obtained with quenched calculations at a comparable pion mass [9] , 
which find an infinite- volume extrapolated value of 5^2 ( 4 He) = 27.7(7.8) (5.5) MeV, close 
to the experimental value. 



B. 1= §, r = 0+ : \Re and 

In nature, the ^He hypernucleus has been well studied experimentally and theoretically. 
The A-separation energy of the ^He J 71 = + ground state is measured to be S\ = 
2.39(0.03) MeV, and for the J n = 1+ first excited state is S A = 1.24(0.05) MeV. These 
two lowest-lying states are consistent with the A coupled to a 3 He J w = ^ + core. A recent 
review of this system can be found in Ref. [51]. 

We have calculated correlation functions in the J n = + channel, which should provide 
the ground state, but not the nearby J n = 1 + first excited state. The sources employed to 
produce the correlation functions are elements of the same 28 irrep of SU(3) as those of 4 He, 
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and hence the extracted states have the same energy 10 . The EMPs from these correlation 
functions are the same as those shown in fig. 14, from which the energies of the lowest lying 
states have been determined, and are the same as those in Table XV. The spectrum in 
this channel, and a subset of associated continuum states, are the same as those in fig. 15. 
There are no continuum states from other SU(3) irreps lying lower than those associated 
with the 4 He spectrum (assuming that we have correctly identified the ground states in the 
three-body sector). However, due to the presence of different SU(3) irreps in this channel, 
the spectrum of excited states of the nucleus, and the continuum states, is expected to be 
different from that in the 4 He channel. 

As is the case for 4 He, while the lowest-lying state extracted from the correlation functions 
has a central value that is lower than any of the non-interacting continuum states, the 
precision of the calculation is not sufficient to completely exclude the possibility that it is 
a continuum state, e.g. 3 He+A, or A He+N. The extrapolated binding energy is given in 
eq. (15). 

C. I = 1, J n = 0+ : AA He, AA H, and nnAA 

At the SU(3) symmetric point, with a deeply bound H-dibaryon, bound di-neutron and 
attractive An interaction, we naively expect to find that AA He and its isospin partners are 
bound. This is in contrast to the situation at the physical point, where a doubly strange 
hypernucleus that is stable against strong decay has not been conclusively observed (for 
recent reviews of the status of experimental investigations into doubly strange hypernuclei 
see, for example, Ref. [51-53]). The states in AA He (with s = —2 and / = 1) and its isospin 
partners can reside in the 27, 28, 35, 81, 64 and 125 irreps of SU(3). However, the sources 
employed in this work produce correlation functions in the 28 and 27 irreps only, and 
therefore the complete spectrum cannot be definitively determined. EMPs from one of the 
correlation functions are shown in fig. 16, from which the energies of the lowest-lying states 
have been determined. The extracted ground-state energies, only calculated for the system 
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FIG. 16: The EMPs associated with one of the eight J n = + AA He correlation functions computed 
with the 24 3 x 48 (left), 32 3 x 48 (center) and 48 3 x 64 (right) ensembles, with momentum 
|P| = 0. The inner (darker) shaded region corresponds to the statistical uncertainty of the extracted 
energy, while the outer (lighter) shaded region corresponds to the statistical and fitting systematic 
uncertainties combined in quadrature. 



The s = —1, I = h systems of various spin configurations have components transforming in the 81 and 
125 irreps that are inaccessible to our operator construction, but that may in principle contain the ground 
state of this system. 
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TABLE XVI: The calculated binding energies in ^He. "g.s." denotes the ground state. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. 



AAHe 


24 3 x 48 


32 3 x 48 


48 3 x 64 


g.s. (MeV) 


157(7)(22)(2) 


154(14) (19) (2) 


156(16)(21)(2) 



at rest in the lattice volume, are given in Table XVI, and shown in fig. 17. The energy of the 
lowest state in the correlation function with contributions from the 28 and 27 is found to 
be the same within uncertainties with that from a pure 27 correlation function. The energy 
of the lowest state in the 28 is that of the ground state of 4 He by SU(3) symmetry, and 
is significantly larger than that of the 27, and clearly the 27 is dominating the large-time 
behavior of the mixed correlation function. Using the result obtained on the 48 3 x 64 
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FIG. 17: The bound-state energy levels in the J w = + y^He (/^\H and nnAA) sector. The 
points and their associated uncertainties correspond to the energies of the states extracted from 
the correlation functions with the quantum numbers of the ground state of The excited 

state of the ^He , in the 28, has the same energy as the ground state of 4 He. The locations of 
the energy-levels associated with non-interacting A- A He, NA-NA, H-dibaryon-di-nucleon, NA-N-A, 
di-nucleon-A-A, H-dibaryon-N-N, and A-A-N-N continuum states, determined from the two-body 
binding energies given in Table VII and the three-body energies given in eq. (9) and eq. (12), are 
shown. 
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TABLE XVII: The calculated binding energies in A ~ 5 H. "g.s." denotes the ground state. The first 
uncertainty is statistical, the second is the fitting systematic and the third is due to the lattice 
spacing. 



4 He 


24 3 x 48 


32 3 x 48 


48 3 x 64 


g.s. (MeV) 


273(19) (39) (3) 


255(25)(37)(3) 


245(28)(81) 



ensemble as an estimate of the binding energy in infinite volume, we find that 

5 (oo) ( AA He) = 156(16)(21)(2) MeV . (16) 

The ground state is more bound than any continuum state (although we have been unable 
to cleanly isolate the ground state of the doubly strange three-body hypernuclei) and we 
identify this as the ground state of the AA.He, aa.H, nnAA isotriplet. However, it is possible 
that this is an excited state of the nucleus, with irreps other than the 28 and 27 containing a 
lower-energy state. Further, it is also possible that this state is a continuum scattering state 
associated with N+ AA H. Clearly, further calculations are also required to unambiguously 
distinguish the energy of the 27 ground state from that of the 28 excited state. 

VII. FIVE-BODY SYSTEMS 

There are a plethora of five-body systems that can be explored theoretically at the SU(3) 
symmetric point, dictated, in part, by the product of five 8's, 

8 <g> 8 <g> 8 <g> 8 © 8 = 32 1 © 145 8 © 100 10 © 100 TO © 180 27 © 20 28 © 20 28 

© 100 35 © 100 35 © 94 64 © 5 80 © 5 80 © 36 81 © 36 81 
©20 125 ©4 154 ©4 154 ©216 . (17) 

In this work, we explore one five-body state that can be produced by local quark-level 
operators, involving only their upper components, with all five baryons in a relative s-wave. 
Unfortunately, this system, with s = —3, has not been experimentally observed. 

A. I = 0, J T = | + : AZ°pnn 

The KEPpnn state has I — 0, s — —3, J 71 = 3/2 + , and belongs to a 10 irrep of SU(3). 
Extending the standard hypernuclear nomenclature, it may be referred to as A= 5 H. Ex- 
perimentally, it is not clear how such a state could be produced and, given the two-body 
interactions, it is not expected to be bound at the physical values of the light-quark masses. 
The EMPs for this system in each of the lattice volumes are shown in fig. 18, from which it 
is clear that the lowest state is negatively shifted with the energies given in Table XVII. 
It is not clear that the 10 contains the ground state of the system, or if it corresponds to a 
continuum state. 

While it is interesting to study this state for algorithmic reasons, the states of more im- 
portance are those that can be accessed experimentally, those with s = 0, —1, —2. These 
more interesting systems have baryons in a relative p-wave, i.e. p-shell nuclei and hypernu- 
clei, and require retaining the lower components of the quark fields in the local operators 
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FIG. 18: The EMPs of the single correlation function for the AE°pnn state. The inner (darker) 
shaded region corresponds to the statistical uncertainty of the extracted energy, while the outer 
(lighter) shaded region corresponds to the statistical and fitting systematic uncertainties combined 
in quadrature. 



by parity considerations. Unfortunately, we find that such operators have poor overlap onto 
such systems, and produce noisy correlation functions. These nuclei can be accessed with 
nonlocal operators and are the subject of future work. 



VIII. SUMMARY AND CONCLUSIONS 

We have presented the results of Lattice QCD calculations of various of the lightest nuclei 
and hypernuclei with A < 5 and with light-quark masses at the (unphysical) SU(3)-flavor 
symmetric point equal to the physical strange quark mass. These calculations were per- 
formed in three lattice volumes with spatial extent 3.4 fm, 4.5 fm and 6.7 fm, and with one 
lattice spacing of b ~ 0.145 fm. Using a new algorithm to perform the Wick contractions, 
ground-state energies of a number of nuclear states were determined from one or more cor- 
relation function(s) generated from local quark-level operators for systems at rest or moving 
in the lattice volumes. A summary of the binding energies determined in this work can be 
found in Table XVIII, and is shown in fig. 19. The approximate binding energy per baryon, 
which is seen to be significantly larger than found in nature, is also shown in Table XVIII. 

In contrast to QCD with the light-quark masses at their physical values, at the SU(3) 
symmetric point all two-body channels except possibly VS( 3 S'i) contain a bound state in 
their spectrum. The SU(3) 1 H-dibaryon is the most deeply bound two-body state, and its 
excitation, transforming as a 27 of SU(3), is also bound. The nature of the sources used in 
this work, each derived from the same light-quark propagator, are such that states in the 
symmetric 85 of SU(3) are not produced in the correlation functions, and as such, we are 
unable to locate these states in the two-body spectrum. The energy splitting between the 
deuteron and the di-nucleon is found to be smaller than the splittings to the other SU(3) 
irreps, consistent with what is found in nature, and the result of a large- N c analysis. It 
is interesting to note that the deuteron remains a finely tuned system even at this heavy 
pion mass. In nature, the ratio of the deuteron binding momentum to the pion mass (which 
defines the range of the nuclear force) is ^/ M^B^/m^ ~ 0.33, where is the nucleon 
mass and is the deuteron binding energy. This quantity is exploited as an expansion 
parameter in the low-energy effective field theory description of nuclear interactions [54]. 
Our calculations reveal that ^M^B^jvn^ ~ 0.24 at m n ~ 800 MeV, which, by this measure, 
is even more finely tuned than at the physical light-quark masses. 
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TABLE XVIII: Summary of the extracted ground-state binding energies of the nuclei and hyper- 



nuclei studied in this work. 



State 


A s I J n SU(3) irrep Binding Energy (MeV) ~ B/A (MeV) 


d (deuteron) 
nn (di-neutron) 
n£ 

H (H-dibaryon) 
nS 


2 1+ TO 19.5(3.6) (3.1) (0.2) 10 
2 1 0+ 27 15.9(2.7) (2.7) (0.2) 8 
2 -1 | 1+ 10 5.5(3.4) (3.7) (0.0) 3 
2 -2 0+ 1 74.6(3.3) (3.3) (0.8) 37 
2-2 1+ 8^ 37.7(3.0) (2.7) (0.4) 19 


3 He, 3 H 
3 v H(hypertriton) 
3 v H(hypertriton) 
A He, ^H, nnA 
|He 


3 | \ + 35 53.9(7.1)(8.0)(0.6) 18 
3 -1 \ + 35 53.9(7.1)(8.0)(0.6) 18 
3 -1 | + TO 82(8)(12)(1) 27 
3 -1 1 | + 27 69(5)(12)(0) 23 
3 -1 1 | + 27 55(6)(10)(1) 18 


4 He 

l He > A H 

AA He > AA H - nnAA 


4 0+ 28 107(12)(21)(1) 27 
4 0+ 28 107(12)(21)(1) 27 
4 0+ 27 156(16) (21) (2) 39 
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d nn 3 He 4 He nE |He |He ^He H-dib nH AA 4 He 

FIG. 19: A compilation of the nuclear energy levels, with spin and parity J n , determined in this 
work. 

In the three-body sector, we are able to cleanly identify the J 71 " = | + ground state 
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of 3 He and its isospin partner 3 H, and the total binding energy is determined to be 
53.9(7.1) (8.0) (0.6) MeV. In the case of the hypertriton, \R, the states in both the .r = \ + 
and J n = | + channels are consistent with being bound nuclear states and not continuum 
states. They are both found to be deeply bound, with the J n = | state being somewhat 
more bound than the J n = | + state. This is in contrast to the situation in nature, where 
the J 71 " = | + hypertriton is found to be very weakly bound. The J 71 = ^ + ground state of 
|He, and its isospin partners and nnA, are cleanly identified, with a binding energy of 
69(5) (12) (0) MeV, which is substantially lower than the corresponding continuum states. 
Further, the J n = | + |He ground state is observed to be more bound than continuum 
states but is somewhat less phenomenologically interesting, as it does not contain an NNA 
component. 

In the case of 4 He, a bound J 71 " = + ground state has been identified, which, while lower 
in energy than any of the continuum states, cannot be unambiguously identified as a bound 
4 He nucleus because of the precision of the calculations. As the sources employed for ^He 
and 4 He are in the same SU(3) irrep, their spectra are identical in the present calculations, 
and as such, this ambiguity is present for ^He also. The ground state of \ A B.e and its isospin 
partners A \R and nnAA can be clearly identified, with a binding energy of 156(16) (21) (2) 
MeV. 

Finally, we have calculated correlation functions in an exotic five-baryon channel, with 
s = —3. Significantly more calculations will need to be performed in order to cleanly identify 
a ground state in this system, but this calculation has demonstrated that the contractions 
for five-body systems can now be performed. 

It is now clear, but hardly a surprise, that the spectrum of nuclei and hypernuclei change 
dramatically from light-quark masses at the SU(3) symmetric point to the physical point. 
While we had already learned this from the recent work on the H-dibaryon, and nucleon- 
nucleon scattering lengths, this has now been demonstrated to be true for even larger sys- 
tems. While the binding energy per nucleon of the deuteron (and di-neutron) is about 
10 MeV, for 3 He and 4 He it is near 25 MeV. These values are significantly larger than 
the 1.1 MeV, 2.6 MeV and 7.0 MeV, respectively, at the physical pion mass. It will 
be interesting to learn how the various thresholds for binding evolve with the light-quark 
masses. Providing accurate binding energies for any given light-quark masses will require the 
inclusion of electromagnetic effects, the leading contributions of which can be determined 
at the classical level and simply added to the results of the LQCD calculations. A deeper 
understanding of the origin of the binding energies calculated in this work will require a 
series of nuclear few-body calculations, which are beyond the scope of the present work. 
In particular, it is important to understand the relative contribution from the two-body, 
three-body, and higher-body contributions to the A > 3 nuclei and hypernuclei, which can 
only be accomplished using modern few-body techniques. 

Our results suggest that quenching in LQCD calculations produces significantly larger 
errors in the binding of nuclei than it does in the hadron masses. This is not too surprising 
given the modifications to the long-range component of the nucleon-nucleon interaction 
due to quenching. It was shown in Ref. [55] that the hairpin interactions that arise in 
quenched and partially-quenched theories generate exponential contributions to the nucleon- 
nucleon interaction in addition to the usual Yukawa interactions at long distances. Therefore, 
one anticipates significant modifications to the binding of nuclei, especially for finely-tuned 
systems. 
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By diversifying and refining the source structure used to generate the correlation func- 
tions, the continuum states in each channel can be explored. In the case of two-body 
continuum states, such as n+ 3 He in the 4 He channel, the established scattering formalism 
of Liischer will allow for the scattering phase shifts in n+ 3 He to be rigorously determined 
from QCD below the inelastic threshold. For the three-body and higher-body continuum 
states, further formal developments are required in order to rigorously determine multibody 
S-matrix elements. 

Lattice QCD has evolved to the point where first-principles calculations of light nuclei 
are now possible, as demonstrated by the calculations at unphysically heavy light-quark 
masses presented in this work. The experimental program in hypernuclear physics, and 
the difficulties encountered in accurately determining rates for low energy nuclear reactions, 
warrant continued effort in, and development of, the application of LQCD to nuclear physics. 
Clearly, calculations at smaller lattice spacings at the SU(3) symmetric point are required in 
order to remove the systematic uncertainties in the nuclear binding energies at these quark 
masses. While not providing quantities that can be directly compared with experiment, these 
calculations provide valuable information about the quark-mass dependence of spectrum of 
the lightest nuclei, and hence the nuclear forces, and will shed light on the fine-tunings that 
are present in nuclear physics. In order to impact directly the experimental program in 
nuclear and hypernuclear physics, analogous calculations must be performed at lighter quark 
masses, ideally at their physical values. 
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Appendix A: Casimirs of SU(3) 

In order to classify the states of the nuclei into irreps of flavor-SU(3), the quark- level sources 
that generate the nuclear correlation functions are acted on with the quadratic and cubic 
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TABLE XIX: The values of the quadratic and cubic Casimir operators in SU(3), c?,{m,n) and 
c 3 (m,n). 



irrep 


m 


n 


C2 


c 3 


1 














3 
3 
6 
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4 
3 
4 
3 
10 
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35 
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12 


18 


35 
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4 


12 


-18 


64 


3 


3 


15 





81 


5 


2 


20 


30 


81 
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20 


-30 


125 


4 


4 


24 






Casimir operators of SU(3), 

C 2 = E T a f a , C7 3 = E d «bc f a f b f c . (Al) 

a abc 

The Casimir operators acting on an irrep of SU(3) that has a tensor representation with m 
upper and n lower indices, Q^'^Z °^ dimensionality 

d(m,n) = hm + l)(n + l)(m + n + 2) (A2) 

have eigenvalues 

C2(m,n) = -(m 2 + n 2 + mn) + m + n 

c 3 {m,n) = — (2m + n + 3){2n + m + 3){m - n) , (A3) 
18 

the values of which are given in Table XIX for the relevant irreps. 

Appendix B: The Expected Continuum States in the Finite Lattice Volumes 

Given the single-hadron and two-body energies that have been extracted in Sees. Ill and 
IV, the continuum states that are expected to arise in the three-body sectors with given 
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quantum numbers can be estimated. Similarly, the information obtained for the three-body 
systems extracted in Sec. V allows for an estimate of the continuum states in the four-body 
sector, and so forth in higher-body systems. In the figures in the main text, this information 
has been presented as the infinite-volume thresholds for the various possible continuum 
channels. Here, we present an example of the expected spectrum of states in the 4 He system 
in the different lattice volumes used in this work. 

For a non-interacting two-component system, comprised of nuclei A\ and A 2 , the indi- 
vidual components have only back-to-back momenta, 

e { Zm = \/ M I + IpI 2 + ^J M l + IpI 2 • (bi) 

For three or more cluster continuum states (for example d+p+n in the 4 He channel), labeling 
the clusters A 1 , A 2 , . . . , A n , the system has energies permitted by momentum conservation 

E£Z.,A n = ^ it P*) t V^l + IPiP > (B2) 

with the obvious generalization to systems with a nonzero center-of-mass momentum. These 
considerations ignore the interactions between the clusters, which will modify the position 
of the corresponding energy levels. For two-body clusters, it is expected that there will 
be 0(1/ L 3 ) shifts in the continuum energies, but for higher-body clusters the form of the 
energy shifts is not known. In fig. 20 we present the expected (ignoring interactions) FV 
energy levels in the 4 He sector for each of the volumes used in this work. 

With more accurate LQCD calculations and additional interpolating operators, we aim to 
investigate these states in the future. However, this makes clear the difficulty in extracting 
excited states in nuclei from this type of calculation. The continuum states rapidly accu- 
mulate as the lattice volume becomes large, and isolating nuclear excited states above the 
lowest-lying continuum states will be challenging with current technology and algorithms. 
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FIG. 20: Expected energy levels in the J w = + 4 He sector. The blue, green and red lines in 
each column denote the location of non-interacting continuum levels in the 24 3 x 48 , 32 3 x 48 and 
48 3 x 64 ensembles, respectively. The location of the states in the 24 3 x 48 and 32 3 x 48 ensembles 
have been displaced slightly for demonstrative purposes. 
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